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Abstract 

The semi-analytical 0{ag) expression for the renormalization group (5 function in the V scheme 
is obtained in the case of the SU{Nc) gauge group. In the process of calculations we use the 
existing information about the three-loop perturbative approximation for the QCD static potential, 
evaluated in the MS scheme. The comparison of the numerical values of the third and fourth 
coefficients for the QCD RG /3 functions in the gauge-independent V and MS schemes and in 
minimal momentum scheme in the the Landau gauge is presented. The phenomenologically oriented 
comparisons for the coefficients of 0(Q:f) expression for the e+e“-annihilation R-ratio in these 
schemes are presented. It is shown that taking into account these QCD contributions is of vital 
importance and lead to a drastic decrease of the scheme-dependence ambiguities of the fourth- 
order perturbative QCD approximations for the e'’'e“ annihilation R-ratio for the number of active 
flavours, nj = 5 in particular. We demonstrate that in the case of QED with A-types of leptons 
the coefficients of the function are closely related to the ones of the Gell-Mann-Low T function 
and emphasise that they start to differ from each other at the fourth order due to the appearance 
of the extra A^-contribution in the V scheme. The source of this extra correction is clarihed. The 
general all-order QED relations between the coefficients of the /3^ and T functions are discussed. 
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I. INTRODUCTION 


The renormalization group (RG) (3 function is one of the basic quantities of the RG 
method, which was developed in the classical works of Refs. 0. i. i- It dehnes the energy 
behaviour of the renormalized coupling constants of the renormalized quantum field models. 
It is known that in the case when the quantum held model under study has the single 
coupling constant, the perturbation theory (PT) expressions for its RG /9-functions depend 
on the choice of the scheme of subtracting ultraviolet (UV) divergences. 

In QED the hrst and the second coefficients of the {3 function are scheme independent 
and were obtained in Ref. [sl from the analytical calculations of the two-loop approximation 
for the renormalized photon propagator performed in Ref. [^. 

In the momentum (MOM) scheme, dehned by subtractions of the UV divergences of 
the photon vacuum polarization function at the non-zero Euclidean point A^, the QED RG 
(3 function coincides with the Gell-Man-Low function T(aMOM); where the expression for 
coincides with the QED invariant charge, uniquely defined by the combinations 
of the Green functions [^. The expressions for the coefficients of the PT series for T(aMOM) 
depend on the number of leptons N. 

For N=l, i.e. in the case of consideration of the electron only, the three-loop term of the 
T function was calculated analytically in Ref. . This result was generalized to the case of 
the arbitrary number N of massless leptons in Ref. . The V-dependent expressions for the 
four- and hve-loop corrections to the Gell-Man-Low function were evaluated symbolically 
in Refs. and respectively. At iV=l the result of Ref. coincides with the similar 
expression, obtained in Ref. [1^ . This feature should be considered as the strong argument 
in favour of the consistency of the complicated analytical hve-loop calculations, performed 
in Ref. (^. 

Another important scheme, which is used in QED, is the on-shell (OS) scheme. In this 
scheme the photon vacuum polarization function is dehned by subtracting UV divergences at 
zero transferred momentum, while the renormalized on-shell masses of leptons are identihed 
with their experimentally measured values. 

In the physical OS scheme the calculations of the (3{c(os) were performed at the three- 
loop level in the work of 0 . The analytical expression for the corresponding four-loop 

3 . ve-loop contribution 

It is in agreement with the result of the work [ij, where this 
term was obtained at V=1 with the help of the concrete RG-relations. The agreement with 
the outcome of the direct hve-loop calculations of Ref. [3 gives extra conhdence in the 


correction was obtained in Ref. 
was obtained in Ref. 


13 


correctness and self-consistency of the results of the complicated computer calculations used 
in Ref. 


10 


The third class of schemes, which we are interested in, is introduced when the dimensional 
is used. These schemes include the minimal subtractions (MS) scheme 


14 


regularization 

[3 and its modihed variants, namely the MS-scheme [3 and the G scheme [3- If is 
possible to prove that for all these modihcations of the MS scheme the coefficients for the 
RG /3 functions coincide in all orders of PT. 


At N=1 the three-loop correction to the /9(ayfg) function was evaluated in Refs. [17|] and 
3 independently (this result had been also presented in the review of Ref. fl^l). In the 
case of the arbitrary N the three-loop contribution to /9(a|^) was obtained analytically in 
Ref. [^. The computation of the four-loop term was completed in Ref. jsj. The five-loop 
correction to the QED (3 function in the MS-like schemes was calculated in Ref. j^. At 
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A^=l this expression coincides with the result of non-direct analysis, performed in Ref. 

It is known that in QCD the MS-like schemes maintain the explicit gauge independence 
of various RG quantities. This property clarihes why in multiloop QCD calculations the 
MS-like schemes are used more often. In QCD the hrst coefficient of the /3-function was 


computed in Refs. [20|,j2l|] and for the number of quarks flavours Uf < 6 turned out to 


be negative. This feature revealed the existence of the asymptotic-freedom property in the 
gauge theory of strong interactions. The two- loop correction to the QCD (3 function in the 
MS-like scheme was analytically evaluated in 


22 , 23 


and is also negative h 

At the three-loop level the QCD jS function was ana^ically calculated in the MS scheme 


in Ref. 26 . This result was confirmed later in Ref. 2 


The four-loop term of th e Q CD (3 
1 and confirmed in Ref. ^9|. For 


function in the MS-like schemes was evaluated in Ref. [_ 

nf=6 the three-loop correction to the (3 function in the MS-like schemes is positive (see the 
numerical results presented below). Note however that the resummation of the PT series 
for the QCD (3 function in the MS-like schemes gives the arguments that this feature does 
not affect the asymptotic freedom property [s^. 

In QCD one can also use anoth er g auge-independent scheme, namely the V scheme. It 
was first introduced in Refs. j3l[, |32| and is determined by perturbative high order QCD 
corrections to the static potential. This scheme was used in Ref. [33| to model massive 
dependence of the first two coefficients of the RG /? function in the V scheme and for the 
related analysis of the manifestation of the massive-dependent corrections in the effect of 
running of the QCD coupling constant from the energies above the production of charm 
quarks to the high energy region above the scale s=M|. In this case the advantage of 
using the V scheme and not the MS scheme is contained in the possibility of modelling the 
smooth transition of the QCD coupling constant through the thresholds of heavy quarks 
productions. Among other applications of the V scheme in QCD is the analysis of the 
perturbative QCD predictions for —)■ bb) [s^. It was shown in this work that within 

the large /3o-expansion the perturbative approximations for r{H^ —)■ bb) in the V scheme are 
converging to the concrete stable value faster than in the MS scheme. 

However, to analyse more carefully the behaviour of various perturbative QCD series for 
the observable physical quantities in the V scheme it is necessary to know high-order PT 
corrections to the QCD f3 function in this scheme. In the present work we will get the 
semi-analytical result for the fourth coefficient of the QCD [3 function in the V scheme, i.e. 
for the (3^ function. In Sec. II the available results of the analytical and semi-analytical 
calculations of the PT QCD corrections to the static potential VQCD(«s(h^)) in fhe MS 
scheme are summarized. The concrete three-loop results, obtained by two groups of authors, 
are compared. Section III is devoted to the definition of the V scheme and to the presentation 
of the concrete results for the third and fourth coefficients of the /3^ function. The problem 
of finding the analytical expression for the concrete known numerical contributions to the 
fourth-order term of the /3^ function is raised. In Sec. IV the numerical values for the 
scheme-dependent coefficients of the QCD (3 function in the V scheme are compared with 
the similar terms, obtained in the MS-like schemes and in the gauge-dependent minimal- 


MOM (mMOM) scheme widely used at present, defined in Ref. [3^. We also get the 0{a 


^ Its first calculation [2^ contained a bug, which resulted in the positive value of the two-loop term and 
in the appearance of the IR-fixed point of the two-loop PT approximation of the QCD fj function. This 


unexpected conclusion stimulated recalculations of this scheme-independent correction 22|- 2^ . They 
resulted in the disappearance of the perturbartive scheme-independent IR-hxed point in QCD. 


3 




















expression for the e’''e“-annihilation R-ratio in the V scheme and compare it with gauge- 
independent scheme and gauge-dependent mMOM scheme results. In fact both MS and 
mMOM scheme were applied recently for the analysis of the behaviour of R-ratio in the 
fourth order of PT Ref. 3^|. Using the concrete physical input, we modify this analysis and 
emphasize that it is more consistent to perform this comparison for Uf = 4,5 numbers of 
active flavours in the energy-region above the production of charm-quark pairs and below 
s ~ 900 GeV^, where the effects of subprocess —?• —)■ hadrons did not yet start to 

manifest itself. In Sec. V we consider the QED limit of the results obtained in Sec. Ill and 
obtain the expression for the 0(0^) approximation of the function in QED. The origin 
of difference with the QED Gell-Man-Low T function, which is starting to manifest itself 
from the fourth term, is demonstrated and explained. The existing common features of the 
PT series for the T and functions is clarihed in all orders of PT. 


II. PRELIMINARIES: THE HIGH-ORDER EXPRESSION OF THE STATIC PO¬ 
TENTIAL IN QCD IN THE MS SCHEME 


Let us first summarise the available information about the perturbative QGD contribu¬ 
tions to the static potential known at present. This physical quantity is used in various 
phenomenologically oriented QGD studies, e.g. in the process of theoretical determinations 
of the charm, bottom and top quark masses, and for the studies of the properties of different 
mesons, composed from the c and b quarks (see e.g. (37| , |38| , |39| and references therein). 

Within PT the static potential in QGD is defined as a renormalized expression for the 
potential of interaction at a distance r between static heavy quark Qh and anti-quark Qf^. 
It is expressed through the following Fourier representation 

d?q 


VQCD(Ai r , «s(/i )) — 


(27r)3 

d?q 


e"^W(g2,/i ,a,(p )) 


( 2 . 1 ) 


^iqr 


— dvrGi 


Ois,V 






where Oisy{q^/p^v) is the renormalized QGD coupling constant in V scheme, Os/dvr = 
g^/167r^, g is the strong coupling constant of the QGD Lagrangian, T“ is the generator 
of the SU{Nc) group, normalized as T“ = A“/2 and Cp is the Gasimir operator, defined 
as = CpSij. In the V scheme its coupling constant asy{q^/Py) is related to the 

numerator of the momentum representation of the static potential in the MS scheme dehned 
in Eq. (12.Ih and is expressed as 


otsy{q^/Pv) = as{p‘^)P{as{p‘^),L) = a^(/i^) V P,f®(L) 


MS/- T\ f ^s{p‘^) 


71=0 


dvr 


( 2 . 2 ) 


The rhs of Eq. (12.2p is expressed through higher-order PT QGD corrections to the static 
potential P^^{L) in the MS scheme which are known at present up to 0(af)-level and will 
be presented below. 

The evolution of the MS scheme coupling constant as{p^) (which depends on the MS 
scheme renormalization parameter p^) is governed by the QGD MS scheme (3 function: 


P 


2 


djas/^Ti) 

dp‘^ 



OO 





i+2 


(2.3) 
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where = asj^Ti and its known fonr MS scheme coefficients, taken from the work of Ref. 
28|, read: 


1^0 — ~ 3^^^^ 

/9i = y C-l - ^CpTpni - y 


j3^ = , 0^2 


-Cl + 2Ci^TFni 


205 

“ 9 “ 


54 

rp2 2 , ^58 ^ 2 2 

~^CFTpTii H XTpTii 


CfCaTfUi 


1415 

~W 


C^TfUi + 


9 

/150653 44 

“ I 486 y 


27 

C(3))C1 + 


39143 136,,„, - . 

y ^ y“C(3) ) CjFpni 


/7073 656 > 2 ^ ^ ^ 

V W “ “9“^^^^ > C^CfTfUi + 


4204 352 \ F^2rp 

H :^C(3) 1 CaCpTfTLi 


27 


9 


AGCfTfUi + 


/7930 224,, X \ ^ 2^2 2 

‘ + —C(3) ClT^Fnl 


V 81 


9 


/1352 

lyr 


^C(3)) Cinn] + ( 


/■ 17152 


243 


9 
704 


424,^ . o 1232,^ . , / 80 

- CaTIu^ H- CfTIu} +-^-C(3) 

243 ^ ' 243 ^ ^ ^ I 9 3^^^ 


^abcd^abcd 


512 1664^^^^! 

“9 ITa 


rii + 


™ + 5tHc(3) 

9 3 ^ 


Na 

^abcd^abcd 


Na 


-n. 


(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 


The characteristic colonr strnctnres of the gronp SU(Nc) are defined as in the detailed work 


of Ref. j40|]. In the notations of Ref. (40| we have [T°’,T°] = where are 


the antisymmetric (nnder permntations of any pair of indices) strnctnre constants, which 
satisfy the well-known relation y“cdy6cd _ Casimir operators, 

Tr(T“T^) = Tf^'^^, Na is the nnmber of the generators of the Lie algebra of the SU{Nc), ni 
is the nnmber of qnarks flavors, d^F^ = Tr{T°'T^^T'^F^)/Q is the totally symmetric tensor. 
The notations are dehning the procednre of symmetrisation of the generators 
^abcd _ Tr{C°‘C^^C^C^'>)/Q is the total symmetric tensor of {C‘^)bc = —where (7“ are 
the generators of the adjoint representation of the Lie algebra of the S'-gronp. The 
corresponding colonr strnctnres in Eqs. f|2.4p - fl2.7p have the following form |28|: 


Ca = N, 

^abcd^abcd 

Na 

^abcd^abcd 

Na 


Cf = 


Nj - 1 


2iV. 


Na =N: -I 


Nl{Nl -t- 36) df^^df^^ N^{Nl -t- 6) 


24 

Nt - QNl + 18 
96iV2 


Na 


48 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 
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The terms, proportional to the n-th powers of L = In (in the polynomial 


of Eq. fl2.2p are expressed as P^^ = 1, P^^{L) = + (3oL, P 


MS 


_ _ _ _ _ _ _ (L) = + 

(2afS/?o + /9l)^ + /5o^^ Pi^HL) = af^+ {3af^^o + ‘2af^^i +(3f^)L + {3af^(3^ + lPo(3i)L^ + 
I3qL^. The powers of L in the expressions presented above arise from the solutions of the 
corresponding RG equations in the MS-like schemes at the three-loop level. 

The coefficients are calculated from the concrete Feynman diagrams. The hrst one, 


,MS 


, was calculated long time ago in Refs, and has the following form 


.,MS 31 


20 . 


1 = — Ca - Tf'iT'I 

1 9 ^ 9 ^ ' 


( 2 . 11 ) 


where ni = Uf — 1. The coefficient was obtained in the 3l|. The bug in the pure 
Yang-Mills contribution to evaluated in Refs. M|, was detected in Ref. 


The final result of these analytical calculations of Refs. [Mj], [32 


IS 


= 


4343 , „ TT^ 22 

162 4 3 

IS \ 

--16((3)]Cj.T^ni + 


C( 3 ) |C 5 - + yC( 3 )) CaTpu, 


20 

~9 


V 81 

2 


Pf^i 


( 2 . 12 ) 


The three-loop constant perturbative contribution to the static potential in the MS scheme 
can be presented as 

^ .. (2.13) 


,MS _ ^(3).„3 , ^( 2)^2 


( 1 ), 


03 =03 nf -|- O 3 Til + O 3 nt + O 3 


,( 0 ) 


The ? 7 ,;-dependent terms were computed in Ref. and have the following form: 
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,(3) 


Oq — 


( 2 ) 

Oo = 


ol^^ = 


20 

~9 


rjid 

-^F 


/12541 368..,^, 647r^ 

+ ^C(3) + 


V 243 

-709.717GlTi. + 


CaT^ + 


/14002 


^C(3) I CfT. 


286 

~9~ 


'A 

296 


135 

71281 \ 

^g2—^ 264((’(3) -|- 80((’(5) ) CaCfTf 


\ Jabcd Jabcd 

g C(3) - 160C(5) j ClTp - 56,83(1)^^^ 


(2.14) 

(2.15) 

(2.16) 


where the error of numerical calculation of the G^Tp-coefficient in Eq. fl2.16p is not indicated 


in Ref. 43 


It is worth emphasizing that in the QED limit with Ga= 0, the analytical expressions of 
the n/-dependent terms, which are proportional to the powers of Tp in Eqs. fl2.lip . (12.12 ) 
and in Eqs. fl2.14p - fl2.16p . are in agreement with the MS-scheme results presented in [44 1 



^ This correction was confirmed later by the author in Ref. 


l3l| 
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The numerical expressions of the nrindependent contributions to Eq. fl2.13p were ob¬ 


tained in Ref. 45 and read 


.( 0 ) 


= 502.24(l)C'l - 136.39(12) 


^abcd^abcd 

Na 


(2.17) 


These results should be compared with the results of the independent calculation of Ref. 


46 


4°^ = 502.22(12)C'l - 136.8(14) 


^abcd^abcd 

N. 


(2.18) 


which have greater inaccuracies. Recently the more accurate result for the second term in 
Eq. fl2.18p was obtained with the help of the computer code used in Ref. 46| . The improved 
result for Eq. fl2.18p is: 


.(0) 


= 502.22(12)Cl - 136.6(2) 


^abcd^abcd 

Na 


(2.19) 


The numerical expression of the coefficient before the second structure in Eq. fl2.19p is in 
agreement with the numerical expression of the same coefficient in Eq. fl2.17p and demon¬ 
strates the reliability of the computer codes, created in the process of calculations, which 


were performed in Ref. and Ref. 


The three-loop nrindependent correction to the static potential also contains the RG non- 


controllable additional term Sn^C^L It is associated with the infrared (IR) divergences, 
which begin to manifest themselves in the the static potential at the three-loop level 
49| . In the effective theory of heavy quarkonium - nonrelativistic QCD- these IR-divergent 


L-terms are cancelled by the concrete UV-divergent contributions (see e.g 

Among the aims of this work is the determination of the four-loop approximation of 
the RG /9 function in the V scheme. This can be done by application of the RG-motivated 
effective charges (EGH) approach, developed in all orders of PT in the works of Refs. 5^, 51 
and independently at the next-to-leading order (NLO) in Ref. (for the concrete NLO 
applications see e.g. the work [s^) The fourth-order approximation of the /3 function in 
the V scheme defines the evolution of asy in the region of intermediate and UV values of 
energy scales. It does not depend on the manifestation of IR physical effects and on the 
RG-uncontrollable L-dependent corrections to the static potential. In view of this we will 
not consider them in our further analysis. 


III. THE FOURTH ORDER APPROXIMATION OF THE QCD /3 FUNCTION IN 
THE V SCHEME 

A. The scale-scheme dependence ambiguities. 

Let us start this section from writing the RG equation for the static potential, which is 
defined in Eq. fl2.ip . In the massless limit, considered in this work, it has the following form 


We are grateful to Y. Sumino for informing us of this new unpublished result of his personal calculations. 
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In QCD the scheme-dependence feature of the PT series for the RG (3 function is the 
more delicate issue than the scheme-dependence problem of the QED RG (3 function dis¬ 
cussed in the Introduction. Indeed, contrary to the QED case, in this realistic theory of 
strong interactions it is impossible to introduce straightforwardly the gauge-invariant ana¬ 
log of the MOM scheme (see e.g. js^, [^,{5^) and thus to construct the invariant charge in 
a unique gauge-invariant manner. In QGD the number of the invariant-type charges of the 
MOM schemes is proportional to 4, namely to the number of vertexes of the Lagrangian (i.e. 
of the gluon-quark-antiquark, gluon-ghost-ghost, three-gluon and the four-gluon vertexes). 
Moreover, the dehnitions of these invariant-type charges depend on different kinematic con¬ 
ditions for hxing the scales of subtractions of UV divergences in the renormalized Green 
functions, which enter these different QGD invariant-type charges. Indeed, hxing the kine¬ 
matics conditions by a different way it is possible to construct a number of MOM schemes, 
i.e. the symmetric MOM scheme j^, the variant of symmetric MOM scheme with one 
external zero momentum ^ and the asymmetric MOM (AMOM) scheme Different 


gauge-dependent MOM schemes were used in the direct calculations of the massless two-loop 
[3^, three-loop jssj, and even four-loop [s^, j^, corrections to the 

QGD f3 function. These analytical calculations revealed the importance of the careful study 
of the dependence on gauge parameter"^. The classical example of the validity of this state¬ 
ment is the discovery that in the AMOM the non-proper choice of the gauge in the two-loop 
PT correction to the QGD (3 function can destroy the asymptotic freedom property of the 
perturbative QGD j^.js^. 

Summarizing the discussions of the gauge ambiguities in the QGD analogs of the invariant 
charges of various MOM schemes, we stress that in these schemes it is impossible to construct 
gauge-invariant analog of the Gell-Man-Low function. In view of this it is important to 
study the expansions of the (3 function in terms of physical coupling constants, which enter 
the effective LO approximations of the RG-invariant physical quantities, e.g. the effective 


coupling constant of the V scheme dehned by the QGD static potential (33 


In all these studies the EGH method, developed in Refs. (50[-j^, was used. To remind 
the basis of this approach consider Erst the system of Eq. fl2.ip and fl2.2p . which defines 
the expansion of the QGD coupling constant of the V scheme through the QGD coupling 
constant in the MS scheme. 

At the hrst step, following the NLO dehnition of the EGH scheme, we dehne the effective 
scale of the V scheme as _ 

/i^ =exp[afV/^o]/i|js (3-1) 

where ap = f Ga - 


g g TpUi and /^o is the hrst coefficient of the QGD (3 function, defined 

in Eq. fl2.3p . At the next step we hx = /iy in Eq. fl2.2p and get the following relation 
between the effective charge of the V scheme and the QGD coupling constant 




~ '^s,Ms(h'v) 


l + a^ 




dvr 


, MS/^“s>Ms(h'v) 


dvr 


(3.2) 




Now it is possible to dehne the EGH /3 function of the static potential, which is the RG 


^ It is worth emphasizing that in the Landau gauge the two-loop expression of the QCD 3 function in the 
number of MOM schemes coincide with the MS scheme results. 
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f5 function in the V scheme 


2 d{asy/^Ti) 


= {as,v) 




i=0 


( (y-sy 
Y 47r 


i+2 


(3.3) 


where asy = asy/47r. The standard RG equation relates function to the /? function in 
the MS-like schemes: 


/9^Kn(a.,Ms(/^v)) =/5 ““Kms(/^v)) 


MS/ 


^ do,sy{cis,Msi^A^)) 


(3.4) 


Consider now the relation between (3 functions, computed in the gauge-invariant UV sub¬ 
traction schemes: 

(3.5) 

dag 

where we use the similar normalization conditions for both /3{as) and /3{ds) function, namely 


2d{ag/4:7i) 

^ a/i 2 




CXD 





(3.6) 


with dg = dg/An. For these normalization conditions the coupling constant of one gauge- 
invariant renormalization scheme dg{^) is related to the coupling constant ag{fi) of another 
gauge invariant renormalization scheme by the following expression: 


) I 1 “t“ (l\ 


ctsify^ 

dvr 


+ 02 


dvr 


+ 03 


dvr 


+ 0(at) 


(3.7) 


Taking into account Eq. fl3.5p . the definitions for (3{dg) in Eq. fl3.6p and the relation of Eq. 
fl3.7p . it is possible to get the following links between the coefficients of the {3 functions in 


two gauge-invariant schemes: 

j3o = I3q (3.8) 

h = (3.9) 

^2 = /d2 — Oi/?i -|- (02 — al)/3o (3.10) 

(33 = 13^ — 2ai(32 + o^/3i -I- (203 — 601 O 2 -I- 4a^)/?o (3-11) 


These formulas reflect the transformation laws of the (3 function from one gauge-invariant 
renormalization scheme to another one. 


B. The V scheme /3 function in QCD : Its 0(a® y)-approximation. 


Consider now the fourth-order approximation of the QCD [3 function in the V scheme. It 
is related to the QCD (3 function of the MS scheme via the Eq. (13.dp . Its gauge-independent 
coefficients can be obtained from Eqs. fl3.8p - fl3.lip . where 


, (3.12) 

Pi = /3p = - dCFTm - ‘^CaTfu, , (3.13) 
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and j3i = fif ■, (3i = /3f^^ with i = 2,3 and aj = for j = 1, 2, 3. Using the concrete resnlts 
for (with i = 0, 1,2,3) and aj = (with j = 1,2,3) from Eq. fl3.10p and fl3.1ip of 
Sec. II, we get the third and fonrth coefficients /Ij and of the QCD (3 fnnction in the V 
scheme: 
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704 512 

+ ^ 


9 


C(3) 


.pHi 


Jabcd^abcd 

2 

■Til 


Na 


22 

— ■ 136.39(12)Ca 
3 


99 jabcdjabcd 

- • 56.83(1)C,4 ^ ^ ni 

3 ^ ^ Na 

o JabcdJabcd o 

I ■ 56.83(1) TfuI + I ■ 136.39(12) 


^abcd^abcd 


Na 


^abcd^abcd 

Na 


TfUi 


The property of the scheme independence of the coefficients /3Y within the gange- 
independent MS-like schemes is the conseqnence of application of the ECH approach to 
the static potential. Indeed, it is possible to show that these coefficients are related to the 
massless gange-independent scheme invariants, introduced in the work of Ref. (for the 
details of derivation see e.g. Ref. 0 ). The analytical expression for Eq. fl3.14p was ob¬ 
tained in Ref. and agrees with the similar one of Ref. |3l| with the C^-term corrected 
later on. 

The result of Eq. fl3.15p is new. Its semi-analytical form is explained by the similar 
representation presented in Sec. II for the coefficients of Eq. fl2.16p . fl2.17p and of Eqs. 
fl2.18p . fl2.19p . obtained in the works 431, 45\ and by the authors of Ref. 4^| respectively. 

Consider now the real QCD case, based on the SU{Nc=3) gauge group of colour. In the 
fundamental representation its group structures are hxed as Ca = 3, Cf = 4/3, Tf = 1/2, 
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Na = 8, dA^'^d'A’^'^ = 135, d^^'^'^d'A^'^ = 15/2 and d'^'^'^d'^'^'^ = 5/12. Converting now the 
SU (iVc)-group expressions presented above for the coefficients of the QCD function into 
the form corresponding to the SU{3) group, we get the well-known results for /3o and /9i 

/3o = 11 - 0.666666nz , (3.16) 

/?i = 102 - 12.66666/1; , (3.17) 

and the following numerical expressions for the third and fourth coefficients of the QCD /3^ 
function : 


/32^ = 4224.181 - 746.0062/1; + 20.87191n2 , (3.18) 

= 43175.06(6.43) - 12951.700(390)/!; + 706.9658(6)/if - 4.87214nf . (3.19) 

The errors of the hrst three terms in Eq. fl3.19p are dehned as the mean square error a = 

\JY 1 \=i where cxj are the numerical errors that arise from the multiplication of the factor 

2 Bn by the computed errors of the corresponding MS-scheme numbers for and , given 
in Eqs. fiTTOll and fiTTTll . 


C. The guess about analytical representation of the numerical terms in the SU{Nc) 
expression for 

It may be inspiring to make a guess on the possible analytical representations of the 
results of numerical calculations of the concrete terms in and Og^^ coefficients. There 
is the general rule that the rate of transcendentality structure is increasing with increasing 
order of PT calculations. 

Following this general rule and considering the terms in the expressions for B 2 
we claim that the numerically evaluated contributions in the expressions for the Og^^ and Og^^ 
coefficients, which enter the expressions for the concrete terms in B 3 , can be decomposed in 
terms of rational and transcendental numbers in the following way: 


709.717 = Ri + R2n^ + Rsti^ + R^CiS) + R57i^C{3) + ReCib) (3.20) 

502.24(1) = -|- R^tt^ -\- i^gTT^ + i7ioC(3) + 7?ii7r^(/(3) -|- 7?i2C(5) (3.21) 

56.83(1) = Rii -\- i?i27r^ + 7?ig7r^ -|- i?i4(/(3) (3.22) 

136.39(12) = Ri3 -\- i^igvr^ -1- i7i77r^ -|- i7ig(/(3) (3.23) 


where Ri are still unknown rational numbers. Note that the rational number is any number 
that can be expressed as the ratio (p/q) of two integers with non-zero q. Thus, some of Ri 
coefficients in Eqs. (I3.20p - (l3.23p may be zero. There are indications that Ru and Riq may 
really be zero. It will be interesting to check this guess by analytical calculations of the 
corresponding complicated Feynman diagrams. 
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IV. THE APPLICATIONS OF THE V SCHEME IN PERTURBATIVE QCD AND 
THE RESULTS OBTAINED IN THE MS SCHEME AND THE MINIMAL MOM 
SCHEME 


General discussions 


In the last few years the interest in studying the perturbative expressions for the QCD 
(3 function in the gauge-independent and gauge-dependent schemes increased. This interest 
was pushed ahead by the considerations of the purity of the conformal windows related to 
the IR hxed points in the expressions for the (3 functions of the strong interactions theories, 


based on the concrete non-Abelian groups with fermions (see e.g. [65l|-|67l|i 


There are also more phenomenologically motivated studies of the behaviour of various PT 
QCD contributions to the RG-invariant quantities, evaluated in the different UV-subtraction 
schemes. The hrst study of the gauge-dependence of the three-loop corrections to the e’''e“- 
annihilation R-ratio was made within the AMOM scheme in Ref. j^. However, this work 
was based on the analysis of the gauge-dependence of the AMOM version of the O(a^) 
contribution to this quantity containing the bugs, evaluated in the MS scheme in Ref. j69| . 
It is worth recalling that this MS-scheme result was corrected in Ref. 0 and conhrmed in 
Ref. 71| and later on in Ref. 72| . In view of this it may be interesting to clarify the status 
of the gauge dependence of the available 0(af) approximation for the -annihilation R- 


ratio in the AMOM-scheme using the 0{ag) corrections, evaluated recently in Refs. 73|, 74 


Quite recently a similar analysis was done at the three-loop level in different gauge- 
dependent MOM schemes, and at the four-loop order in the mMOM scheme, specihed for 
the case of the Landau gauge 0]. This mMOM scheme was formulated in Ref. [0 and 
already used in the theoretical studies of the behaviour of the gauge-dependent QCD (3 
function for different numbers of fermions flavours Uf (see the works of Refs. [0], [0, 0|). 
In this section we will compare the expressions for the coefficients of the RG (3 function in 
the V scheme obtained in Sec. m with the similar mMOM-scheme results. In the next 
section we will use the results of Sec. uni to study the third and fourth order approximations 
of the e’''e“-annihilation R-ratio in the V scheme and compare it with the results obtained 
in the MS scheme and in the mMOM scheme, which were presented in Ref. 


36 


B. The definition of the minimal MOM scheme. 

Let us hrst briehy review how the mMOM scheme is dehned. Using the standard notations 
for the renormalization constants of QCD in an arbitrary linear covariant gauge namely 

A7 = c“ = v^c“, go = Z^g, Aq = (4.1) 

where 'ip, A“, c“ are the quarks, gluons and ghosts helds respectively, g is the constant of the 
strong interaction, A is the gauge parameter, which is included in the Lagrangian QCD as 
(cI^A“)^/2A. We hrst write down the non-renormalized gluon propagator in the momentum 
space: 


^ab 2 


i5, 


ab 


+ le 


-g^^+{l-X) 


-|- is 


(4.2) 
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The form of the QCD Lagrangian dictates how to relate different renormalization constants. 
For example, the renormalization constant of the gluon-ghost-ghost vertex has the following 
form: 


^ccg 


= z, 


7I/2 7 


(4.3) 


The dehnition of the mMOM scheme is based on the consideration of this relation 
Taking into account Eq. fl4.3p one can write down the expression for the QCD coupling 
constant of the mMOM scheme as 


a 


mMOM 


(f^) = 


^mMOM(/^2)(^mMOM(/^2))5 


-a. 


(4.4) 


Following the proposals of Ref. [3^ the renormalization expressions for the gluon and ghost 
propagators are dehned by using the requirements that at = /i^ their residues are equal 
to unity, namely 




a 


mMOM 


{P^))\r'‘-ii? = 1 


(4.5) 


Then the renormalized expression for the gluon propagator, dehned in the Landau gauge 
A=0, will take the following form 





p^p’' \ 

p"^ + ie J 


^(p2^^mMOM(/^2)) 

-|- ie 


(4.6) 


while the expression for the ghost propagator is dehned as 


^ab 


= id, 


ab~ 


p^ -|- ie 


(4.7) 


The most important additional requirements of the mMOM scheme [3^,j^ are the spe¬ 
cial dehnitions of the renormalization constant of the gluon-ghost-ghost vertex and of the 
renormalization constant of the gauge parameter, namely 


^mMOM(^mMOM) ^ ^MS(^MS) 


^mMOM(^mMOM) ^ ^MSj-^MS^ 


(4.8) 


Taking into account the dehnition of the QCD coupling constant in the MS scheme through 
the same vertex 


a 


MS 


(f^) = 






-a. 


(4.9) 


and Eqs. fl4.3|) and fl4.8p . one can get the useful relations between the renormalization 
constants of the mMOM and MS schemes 


^mMOM ^^mMOM^mMOM _ 


and the relation between the renormalized QCD coupling constants of these schemes 


a 


mMOM 




ymMOM / 7 mMOM\ 2 _ 

'MS 


^MS 


V ZMS ) 


a 




(4.11) 


All formulas written above are valid for any linear covariant gauge and for the Landau 
gauge A=0 in particular. This choice of the gauge leads to the simplihcation of the hnal 
perturbative results we will be interested in. Note also that the application of the Landau 
gauge allows us to simplify dehnite lattice Yang-Mills studies (see e.g. 111 ). 
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C. Comparison of the fourth-order approximations of the QCD /? function in the 
V, mMOM and MS schemes. 



four-loop computations. In the Landau gauge they take the following numerical form 


^mMOM,A-o ^ 3040.482 - 625.3867nz + 19.38330nf (4.12) 

^mMOM,A=o ^ 100541.05 - 24423.330n; + 1625.4022nf - 27.49263nf (4.13) 


It is interesting to compare these results with the numerical expressions of the same coeffi¬ 
cients of the QCD (3 function in the gauge-invariant V scheme (see Eqs. fl3.18p and fl3.19p l 
and in the gauge-invariant MS scheme, namely with 


^ 1428.500 - 279.6111ni + 6.01851nf , (4.14) 

= 29242.96 - 6946.28971; + 405.0890n;2 + 1.49931nf , (4.15) 


26 and 28 


which follow from the results of analytical calculations of Refs. 

For the completeness, in Table 1 we present this comparison for all numbers of quarks 
flavours 1 < tij < 6 , where nf = rz; + 1. These notations are identical to the ones used for 
hxing the numbers of heavy flavours, which are considered in the PT QCD expression for 
the static potential V, where rii is the number of quarks, lighter than Qh- They enter virtual 
corrections among the heavy quark and antiguark of the flavour Uf and vary in the region 
3 < 77; < 5. 


The numerical coefficients of the QCD (3 function in different schemes 




Ay 



nmMOM,A=0 

P2 

nmMOM,A=0 

Ps 

1 

3499.047 

30925.46 ± 6.44 

1154.907 

22703.26 

2434.478 

77715.63 

2 

2815.656 

20060.55 ± 6.48 

893.351 

16982.73 

1867.242 

57976.06 

3 

2174.010 

10551.11 ±6.54 

643.833 

12090.37 

1338.771 

41157.38 

4 

1574.107 

2367.90 ± 6.62 

406.351 

8035.18 

849.068 

27094.64 

5 

1015.948 

-4518.30 ±6.72 

180.907 

4826.15 

398.131 

15622.88 

6 

499.533 

-10136.74 ±6.84 

-32.500 

2472.28 

-14.038 

6577.14 


Table 1. The comparison of the numerical values of the third and fourth coefficients of the 
QCD f3 function in the V , MS and mMOM scheme in the Landau gauge. 


The results of this Table demonstrate that the asymptotic structure of the PT series for 
the effective (3 function in the V scheme has the non-regular behaviour and differs from the 
asymptotic structure of the PT for the (3 function in the MS scheme, which was considered in 


Ref. [76| using the approach developed in Ref. j77|. In view of this it is of interest whether 


this non-regular behaviour of the PT series for the (3^ function will manifest itself in the 
process of studies of scheme dependence of high-order coefficients for the characteristics of 
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typical physical QCD processes, e.g. for the e+e“-annihilation R-ratio in the region of direct 
production of the pair of heavy quarks and antiquarks with nj = 4, 5 numbers of flavours. 
We will not consider in this work the case of n/ = 6 , related to the direct production of the 
pair of tt-quarks in the process e’''e“ —)■ hadrons, which may be studies in future if the ILC 
will be built. Indeed, the total cross section of this process is dominated by the subprocess 
—)■ —)■ hadrons and not by the subprocess e+e“ —hj—h hadrons that interests us in 

this work. 


D. The fourth-order approximation for the e'''e R-ratio in the MS and V schemes 

We now discuss the fourth-order PT expression for the e''“e“-annihilation R-ratio in the 
V scheme. The idea to study this particular expression, as well as the PT expressions for 
other observable physical quantities in the V scheme, was proposed some time ago in Ref. 
78l |. In this section we will realise this proposal, obtain the fourth-order V scheme PT 


approximation for the e’''e“-annihilation ratio -R(s) and compare its coefficients and energy 
dependence with the results, obtained in the MS scheme and in the Landau gauge variant 
of the mMOM scheme 3^. The studies to be made in this subsection supplement the ones 
presented above. Moreover, the results obtained in Sect. IV C will be used in the process 
of the numerical calculations to be presented below. 

We remind the reader that the e'''e“-annihilation R-ratio is dehned as 




(j(e+e —)■ 


^ ^ = i2.im n(s+i£) 


cTo(e+e ^ 7 ^ ) 


(4.16) 


where s is the transferred energy in the Minksowskian region, (To(e'''e“ —)■ 7 ^ = 

47 r^a/( 3 s) is the theoretical normalization factor, 11 (g^) is the QCD expression for the photon 
vacuum polarization function 


n^,2(g ) = (gf,qu - )n(gQ = i d^x e^^^{ 0 \T j^(a;)>( 0 )| 0 ) 


(4.17) 


and = 'Yh Q is the electromagnetic hadronic current. Since the e'''e -annihilation 
/ 

R-ratio is the RG-invariant quantity, it obeys the RG equation without anomalous dimension 
term, namely 




(4.18) 


In the MS scheme the O(a^) approximation for the e'''e R-ratio has the following form 


^Ms _ 3 ^ g2 1 ^ 4 


a 


MS 


47r 




MS I 


MS 


47r 


^MS 

, MS ( 

^ * 47r 


n/MS 

, MS ( 

+ > 47r 


(4.19) 


where the coefficient was evaluated analytically in Ref. [79| and numerically in Ref. [80 


and conhrmed analytically in Ref. The coefficient was analytically evaluated in 

Ref. [^ and conhrmed in Refs. 0 and j^, while the symbolical expression for the non¬ 
singlet and singlet contributions to were obtained analytically only recently in Ref. 
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and ^ 


form as 


respectively. 


The coefficients and can be expressed in the nnmerical 


= -1.84472n/+ 31.7713 , (4.20) 

rp = -0.33139nJ - 76.8085n/ - 424.763 - 26.44355/ , (4.21) 

rp = 5.50812p- 204.1431p +4806.339/1/-40091.67 (4.22) 

+ (49.0568/1/ - 1521.214)5/. 


where the terms, proportional to 5/ = f Qf Q})^ t^e singlet contribntions. 

In the V scheme the PT expression for the e’''e“ R-ratio is dehned as 


= 3 




Q 


1 , a^s,Y , V 

1 + 4—^ + 

47r 


a 


s,V 


47r 


+ P 


47r 




47r 


(4.23) 


/ 

Using the ECH approach of Ref. jlH and the V scheme relations of Eqs. (13.111 . (I3.10|l and 
(13.lip we obtain the following general expressions for rj : 


P = rp-4ap, (4.24) 

rV = rp - 4ap - 2aprV (4.25) 

r3V = rp - 4ap - - (2ap + (4.26) 

and the nnmerical valnes of these coefficients, namely 

rX = 2.59972/1/ - 9.5620 , (4.27) 

P = 0.50749/iJ + 113.6320/1/ - 2054.140 - 26.44355/ , (4.28) 

P = 3.05815/iJ- 144.9455/i5 +3455.279(2)/!/- 20387.90(1.17) (4.29) 

- (39.0881/1/ + 701.466)5/ 


The errors in the valnes of the n / and /i°-terms in Eq. (I4.29p arise from the nnmerical errors 

in the valnes of the /i/ and /i^-dependent constitnents and a® of the coefficient 
dehned in Eq. (l2.16p and Ref. (l2.17p . which enter into the dehnition of throngh Eq. (l4.26p . 


E. The comparison of the fourth order V- , MS- and mMOM-scheme 
approximations for the e'’"e“ R-ratio. 


As the start of the study of the scheme and energy dependence of the e'''e“-annihilation 
R-ratio in different orders of PT in the case of applications of three different schemes we 
hrst present in Table 2 the comparison of the following from Eqs. (I4.27I1 - (I4.29I1 and Eqs. 
(I4.20p - (l4.2ip numerical expressions for three PT coefficients in the V and MS scheme with 
the numerical expressions of the same coefficients, obtained in the Landau-gauge version of 
the mMOM scheme in Ref. fsbl. 
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The numerical coefficients of the R-ratio in different schemes 


'1 

'1 

^mMOM 

'2 

'2 

^mMOM 

'2 

'3 

'3 

^mMOM 

'3 

1 

-6.9622 

29.9265 

-21.9622 

-1966.444 

-528.346 

-1575.567 

-17815.06T1.17 

-36956.12 

-13190.55 

2 

-4.3625 

28.0818 

-19.3625 

-1830.134 

-584.994 

-1467.688 

-14188.58T1.17 

-31536.11 

-8632.68 

3 

-1.7628 

26.2371 

-16.7628 

-1708.676 

-658.171 

-1374.660 

-11244.00T1.17 

-27361.22 

-4748.58 

4 

0.8368 

24.3924 

-14.1631 

-1602.069 

-747.876 

-1296.483 

-9033.31±1.17 

-24310.08 

-1590.24 

5 

3.4366 

22.5477 

-11.5634 

-1475.696 

-819.494 

-1198.540 

-6434.41±1.17 

-20591.03 

1575.00 

6 

6.0363 

20.7029 

-8.9637 

-1369.944 

-913.410 

-1121.218 

-4775.30±1.17 

-18149.16 

3873.49 


Table 2. The comparison of the numerical values of the known coefficients for the 
e^e“-annihilation R-ratio in the V, MS and in the Landau-gauge version of the mMOM 

scheme. 


Note that the values of the coefficients and -with i=2,3 are negative for 

any number of Uf, apart from the case of value at Uf = 5, 6. In the MS scheme 

this feature is related to the manifestation in the expressions for and of the effects 
proportional to tt^, which arise from analytical continuation to the Minkowskian region of 
energies of the PT contributions in the rhs of Eq. (Id.lbh (for a detailed explanation see 
e.g. Ref. (j^). The negative values of the V-scheme coefficients are also related to these 
kinematic effects, but the numerical difference with the negative values of and 
terms is related to the numerical values of the additions contributions to and r^-terms. 
Note that in the case of they are negative (see Eq. (14.2511 ) but in the case of they 
are positive due to interplay among the third huge positive contribution to Eq. fl4.26p and 
other negative contributions to the same equation. Note also that the values of are very 
closed to but this feature does not remain at the fourth order of PT. 

We now plot the energy and scheme dependence of the next-to-leading order (NLO), next- 
to-next-to-leading order (NNLO) and next-to-next-to-next-to-leading order (N^LO) approx¬ 
imations for the function r(s) = R{s)/{3J2f Q}) ~ 1- ^ depends on s = q'^, where s is 
measured in GeV^. The hrst three plots are presented in Fig. 1 for the energy region above 
the threshold of charmonium production and below the threshold of the bottomonium pro¬ 
duction, i. e. in the region where Uf = 4 numbers of active flavours are contributing to the 
expression for r(s). In Fig. 2 the scheme dependence of the NLO, NNLO and N^LO approxi¬ 
mations of the same function are presented in region with nj = 5 numbers of active flavoures. 
More dehnitely, we consider the energy region above the threshold of bottomonium produc¬ 
tion and up to the energies s = 900 GeV^, where the subprocess —)■ —)■ hadrons, 

which starts to dominate near the beginning of the left shoulder of the direct manifestation 
of Z'^-boson in the collisions, can be safely neglected. 

The energy dependence of coupling constant ag = a^®/(47r) of the NLO, NNLO approxi¬ 
mations of the PT expansions of the e+e“-annihilation ratio R{s) in the MS scheme, which is 
presented in Eq. fl4.19p . is dehned through the powers of logarithmic terms L = ln(s/A^^^)- 
as 


NNLO 


1 


I3i ln(L) 
/33L2 

^^^NNLO 


(4.30) 

(4.31) 
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FIG. 1: Scheme dependence of the NLO (left), NNLO (right) and N^LO (bottom) approximations for the 
e+e” characteristic r(g^) = i?(q^)/(3 X)/Q/) ~ 1 the case of n/=4 numbers of active flavours. The dashed 
black curve depicts the variations of the MS approximants. The solid (green) line demonstrate the variations 
of the mMOM-scheme results, while the solid red line shows the V-scheme results. 


where 


Aa 


NNLO 




[/3pn"(L)-/3pn(L)+A/3„-/3j| 


At the fourth N'^LO, hrst studied in Ref. [82|, one has 


^N^LO _ ^NNLO 
Ci „ - (Jj c 


-\- Act 


N^LO 


where the additional correction reads: 


Aa 


N^LO 


a’L* 


I3f I - lii‘’(L) + r 1„^(L) + 2 lii(L) - r ) - ln(L) + /sA 


(4.32) 


(4.33) 


(4.34) 


In the numerical form the expressions for the MS /^-functions coefficients /3i in Eqs. 
fl4.30p - fl4.34p are defined in Eqs. fl3.16p . fl3.17p and Eqs. fl4.14p . fl4.15|) respectively. For the 
concrete numbers oirif flavours their values are given in Table 1. Note that in the analysis of 
Ref. 3^ the same expansion was used for the numbers of active flavours u/=5 and nf=6 and 
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FIG. 2: Scheme dependence of the NLO (left), NNLO (right) and N^LO (bottom) approximations to r{q‘^) 
are presented for nf=5 numbers of active flavours. The variation of the MS-, mMOM- and V-scheme results 
is indicated by the three curves as in Fig.l 


for the value of = 500 MeV, which did not vary from order to order of the MS-scheme 


perturbative expressions considered in Ref. [3^. In the process of obtaining our results, 
presented in Figs. 1 and 2, and keeping in mind physical motivations, discussed above, we 
used n/=4 and nf=5. 

36] the values of the parameters 


Contrary to the studies of Ref. 


'■mMOM 


and 


parameter (that is new to this work) were not hxed, but depend on the choice of both 
rif and the order of approximations. The concrete results for the values of the parameters 
used are presented in Table 3. 
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The numerical values of the Aqcd in different schemes, MeV 


the order of approximation z/ 

MS 

A^/l 

ny 

^Vmom 

4 

2 

350 

500 

625 

4 

3 

335 

475 

600 

4 

4 

330 

470 

590 

5 

2 

250 

340 

435 

5 

3 

245 

335 

430 

5 

4 

240 

330 

420 


Table 3. The dependence of the parameters, used for getting the results of Figs. 1 and 2 
from the n/, z/ (order of approximation), and from the choice of the scheme. 


In the cases of n/=4 numbers of active flavours and z/ = 2, 3,4 the values for given 

in Table 3 are hxed from the results of the hts hts of the Fermilab Tevatron experimental 
data for the xF^ structure function of the neutrino-nucleon deep-inelastic scattering process 
at the of the theoretical FT results, performed in Ref.j^. In the case of rij = 5 

the values of at z/ = 2,3,4 were obtained in Ref. from the related results for 

using the the NLO, NNLO and N^LO matching conditions, evaluated at the NNLO 
in Ref. js^ and 861 and at the N^LO in Ref. j^. The matching point in these conditions 
was hxed by the on-shell b-quark mass values, extracted at different orders of FT from the 
analysis of heavy quarkonium spectrum while taking into account the Fade estimated value 
of the coefficient as from Eq. fl2.13p . obtained in Ref. [87|. These Fade estimates turned 


out to be in satisfactory agreement with the results of direct calculations of the value of as 
obtained later (see Refs. 




we 


may safely use the values for A^ 


(n/=5) 


46|). In view of reliability of the results of Ref. 

from Table 3 for transforming them to the values of 

the scale parameters and in particular. 

In general the scale parameters of the MS , V and mMOM schemes considered in 
Table 3 are related by the following equations: 

aF*" = A<g'"exp|ap(n,)/ft(n,)| , A^m = Ag^<ixp|(rp(n,) - rr™(nt))/4/?o(n,)] . 

(4.35) 

They are derived by means of the ECH approach. We used these expressions to get in 
Table 3 the numerical values of Ay and from the results described above for 

A^\ Combining them with the numerical values for the coefficients /9j(n/) , /3j(n/) and 

in the analogs of Eqs. fl4.3Up . fl4.32p and Eq. fl4.34p . and 
taking into account the expressions for the coefficients r* in r(g^) = Q}) ~ ^ 

in three different schemes, we plot in Figs.l and 2 the energy dependence of r{q‘^) in three 
different orders of FT and three different schemes, namely MS, V and mMOM schemes in 
the case of nj = 4 and n/ = 5 respectively. 
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F. Discussions of the results. 


Considering now the plots of Figs.l and 2 we may conclude that in all cases the PT 
approximants for the function r(s) related to the e'''e“-annihilation R-ratio are converging 
in all schemes. In the MS scheme the rate of convergence of the related PT approximants is 
better than in the V scheme and mMOM scheme. At the NLO the results of the V scheme 
are closer to the mMOM ones than to the results obtained in MS scheme, while at the NNLO 
the situation is reversed - the V-scheme approximations are closer to the MS-ones, while 
the application of the mMOM scheme puts a lower bound on the theoretical expression for 
r(s). However, at the N^LO the lower theoretical bound on the energy dependence of r(s) 
is changed again and the lower bound is now obtained within V scheme. The comparison of 
three approximants for r(s) in the case of consideration of the V-sche me results supports 
the conclusion, made in Sec. IV. C, that the PT approximants in the V scheme have less 
regular behaviour than the MS ones. The results of Table 2 demonstrate the positive feature 
of taking into account 0(a^)-corrections to e+e“ annihilation R-ratio in all three schemes. 
Indeed, the scheme dependence of the expression for the ratio is drastically decreased 
at this level. This is the positive message, which supports the work presented above on the 
inclusion of the 0(0;^) correction in the theoretical approximations in the MS, mMOM and 
V schemes. 


V. THE FOUR-LOOP QED RESULT FOR THE RG /3 FUNCTION IN THE V 
SCHEME 

Consider now the case of QED with N types of identically charged leptons. We will use 
the results of Sec. IIII Bl for the the fourth-order PT approximation of the RG V-scheme 
/9-function of the SU{Nc) colour gauge group theory. Fixing the SU{Nc) group weights in 
Eqs. flTT^ . (Km . flRTTIl and fl3T5|l as Ca = 0, Cf = 1, Tp = 1, = 0, = 1, 

Na = 1 and rif = N, we obtain the following four-loop semi-analytical expression for the 
RG /9 function in QED in the V scheme: 

= 2 Oy -|- AN Oy ^ — 2N -|- (—C(3) + (^-1) 

+ f - 46iV + (104 + ( 3 ) _ 1^C(5) - 5 ■ 56.83(l))Af2 + (128 - ?pC(3))/V=>) + 0{a%) 

where ay = av/47r and N is the number of leptons. Gomparing this result with the four-loop 
approximation of the QED fd function in the MOM scheme, i. e. of the Gell-Man-Low T 
function, namely with 

4 / 64 184 \ 

4/(aMOM) = 2 A^®mom + i ~ (^-2) 

+ f - 46iV + (104 + 5Hf(3) _ + (128 - ^C(3))A'=‘)a|,oM + 0(al,ou) 
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where omom = «MOM/47r, we conclude that in spite of identical agreement at the third order 
of PT the general expressions for the RG QED f3 function in these two different schemes 
are not the same. They start to differ from the fourth order of PT due to contributing 
to the 0{ay) coefficient of the /5^-function of the additional light-by-light-type scattering 
diagrams, which appear in the QED analog of the coefficient in the MS scheme, given in 
Eq. fl2.16p . They enter in the definition of the A^^-term of the coefficient of the V-scheme 
QED /9-function through Eq. (Id.llh . 

It is possible to clarify what kind of N-dependent high-order coefficients of the following 
expression of the QED f3 function in the V scheme 


/^^(av) = 


OO /X i+2 

y ( ay 


i=0 


An 




i=l 1=1 


Q^v \ 
47r / 


i+2 


(5.3) 


will also receive additional contributions and what kind of the N-dependent coefficients of 
the QED /9^ function will coincide with the similar expressions for the T function, which 
we will define as 


( \ 2 OO i y X i+2 

' i=l 1=1 ^ ' 

Using the analogs of Eq. (Id.lOh and (13.111) . which can be derived using the considerations 
of Ref. j^, we arrive at the following relations: 

/Sr"' = + a/37"' (5.5) 

where extra terms in the N-dependent contributions to the coefficients of the QED 

/9^ function appear in the following region of indexes [i, /] = [i > 3, 2 < / < i — 1]. 

In the cases of [i, /] = [i > 3, / = 1 or i] the proportional to coefficients of the and 4/ 
function, defined in Eq. fl5.3p and fl5.4p . are the same. In the case of i=3, which corresponds 
to the totally known for the moment fourth order results, these identical coefficients are 
proportional to N and N^. At the third order the proportional to N-term was analytically 
evaluated in Ref. [M . At the fourth order of PT the proportional to N and terms were 
evaluated in Ref. |8|. For i = A the terms under discussion can be obtained from the results 
of Ref. and read 


= -^ + 128C(3) 

^v[4] ^ ^ ^ ^ 


(5.6) 

(5.7) 


Note that this result from Ref. is in agreement with the multiloop expression for this 
particular contribution to the Gell-Man-Low function, evaluated in Ref. [9l| up to 20 loops 
analytically and numerically up to 100 loops. The scheme-independence of the linear-in-N- 
contribution to Eqs. fl5.3|) and fl5.4p is the consequence of the conformal symmetry property. 


® This observation was made and used in the unpublished work of A.L.Kataev and A.V. Garkusha, see Ref. 
as well. 
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which is valid in QED in the perturbative quenched approximation (for the recent detailed 
study see Ref. 

In the numerical form the scheme-dependent coefficients of the /dq^p-function read: 

= -2iV + 5.19943Ar2 (5.8) 

(3^ = -46iV + 284.818(26)iV2 - 25.42447Ar^ (5.9) 

The analogous expressions for the three- and four-loop coefficients of the QED (3 function 
in the MS scheme follow from the analytical results of Ref. jsf and have the following form 

^ -2iV + 4.88888A^2 (5.10) 

= -46A^ + 82.9753iV2 + 5.06995iV3 (5.11) 

The numerical expressions for the analogous coefficients of the T function (or the QED f3 
function in MOM scheme), which we obtain from the same work of Ref. [^, are 

^2 = -2iV + 5.19943A^^ (5.12) 

Ts = -46iV + 133.2714iV2 - 25.42447Ar^ (5.13) 

Note once more that the first three coefficients of the (3^ function and of the 4/ function are 
the same and start to differ from the fourth order of PT in the following way 

= ^3 _ i5i.54(2)iV2 (5.14) 

This additional contribution arises from the light-by-light-type scattering contribution, 
which is typical to the V scheme. 

For completeness we present the QED expressions for the 0{a^) approximations for the 
4/ and (3q^d functions in the case of N=1 : 

4i(cimom) = l-3333a]y[;Q]yj-|-dajyfQjyi;-|- 3.1994aj^Qj^ — 153.8469aj^Qj^-|-0(a]yjQ]yj) (5.15) 
(^qedM = 1.3333a^ + 4a^ + 3.1994at,- 305.3936(266)a^ + O(a^) (5.16) 

One can observe that even for N=1 the numerical effect of light-light-scattering contribu¬ 
tion, which is typical for the V-scheme (see Eq. f|5.14p ), is rather sizable and almost equals 
to the whole value of the other term in the expression of Eq. fl5.14p . 


VI. CONCLUSION 

In this work we consider the definition of the gauge-independent RG QCD (3 function in 
the V scheme. Using higher-order corrections to the static potential of the quark-antiquark 
interaction and /3 function in MS scheme, we compute the fourth term of the PT expression 
for the (3 function in V scheme in the general case of SU(Nc) group in the semi-analytical 
term. Our guess of possible expressions of the corresponding numerical contributions through 
concrete transcendental numbers is made. The comparison of the numerical expressions of 
the scheme-dependent coefficients of the (3^ function of QCD with the similar coefficients 
of the QCD (3 function in the MS and mMOM scheme in the Landau gauge are presented. 
The indication that the structure of the PT series for the effective (3 function in the V 
scheme has non-regular asymptotic behaviour and differs from the asymptotic PT for the 
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jS function in the MS scheme are presented. The results obtained in the V scheme are 
used to study the scheme dependence of the 0 (q;^) approximation for the e’''e“ annihilation 
R-ratio in the energy region above the thresholds of production of the charmonium states. 
The conclusion is made that the comparison between the fourth-order expressions for the 
e’''e“ annihilation R-ratio, obtained in the MS schemes, in the Landau-gauge variant of the 
mMOM scheme and in the gauge-independent V scheme leads to a drastic decrease of the 
scheme dependence of the fourth-order perturbative QCD predictions for the case of n/ = 5 
numbers of active flavours in particular. Considering the QED limit of the SU(Nc)-group 
function we observe that its perturbative expression is starting to differ from the perturbative 
expression for the Gell-Mann-Low T function from the level of the 0(Q;Y)-corrections. The 
relations between coefficients of the QED function and the T function are presented 
in all orders of PT in the case of the N-types of identical leptons. The conclusion that 
starting from the fourth-order perturbative approximation two N-dependent terms in the 
coefficients of the perturbative expansions of the and T functions will always coincide is 
made. Theoretical reasons of this foundations are presented. 
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